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Abstract 

,-0 \ A survey is given on asymptotic diffusion coefficients of particles in lattices with random 

\ transition rates. Exact and approximate results for single particles are reviewed. A recent exact 

; expression in d = 1 which includes occupation factors is discussed. The utilization of the result 
is demonstrated for the Miller- Abrahams model and a model of random barriers combined with 
random traps. Exact and approximate results for the site-exclusion model in disordered lattices 
are also given. 
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1. Introduction 



Transport processes of particles in disordered materials exhibit a variety of interesting phe- 
nomena, such as a strong reduction of the asymptotic diffusion coefficients, anomalous frequency 
Q\ . dependence of the conductivity, dispersive transport, etc. The explanation of the transport pro- 
cesses of single and many particles in disordered materials has been a challenge to theory; it 
is of great practical importance as well. This paper focuses on the asymptotic diffusion coef- 
ficients of single and many particles for different models of disordered lattices. Quite different 
models of disorder were introduced to describe particle transport in crystals with point defects, 
as well as in amorphous materials and glasses. Exact results will be presented as far as pos- 
sible; these results are mainly but not always restricted to the dimension d — 1. Recently a 
general exact expression for the diffusion coefficient of single particles in d = 1 has been derived 
JT|, 0] . The insight obtained from this derivation can serve as a basis for the effective-medium 
o3 1 approximation in higher dimensions. 

The derivation of diffusion coefficients of many particles where multiple occupancy of sites 
is excluded will be restricted to the coefficient of collective diffusion which appears in Fick's 
law. Exact results for this diffusion coefficient in disordered lattices are scarce. It will be 
shown that a properly formulated effective-medium approximation yields reasonable results in 
higher-dimensional disordered lattices. 

2. Models 

Most theoretical models for the description of transport in disordered materials use ordered 
lattices and put the disorder into the transition rates between the sites. Typically only nearest- 
neighbor transitions are considered. Several models of the disorder have been considered in the 
past. Two prototype models are the random-barrier (RB) model and the random-trap (RT) 
model see Fig. 1 (a) and (b). In the RB model the transition rates between two neighbor 
sites i and j are symmetric, T^- = Tjj. In the RT model, the transition rates Tj originating 
from the sites are given; they are independent of the final sites. 

The transition rates T^, or Tj, respectively, are taken from a common distribution, f(T). 
Usually thermal activation is assumed as the physical process leading to transitions. The sim- 
plest form of a thermally activated process is given by the Arrhenius law, T = Tq exp(—E/kBT). 
In the case of the RB model, the energy E represents the barrier height, while for the RT model 



E is indicative of the trap depth. Instead of specificing the distribution f(T), one may specify 
the energy distribution v(E), from which barrier heights for the RB model, or trap depths of 
the RT model are chosen. Note that all sites have the same energy in the RB model, while all 
barriers lie at the same energy in the RT model. 
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Figure 1: Schematic representation of the four models of disordered transition rates that are 
discussed in the text. 

A third model is the one introduced by Miller and Abrahams (MA) |5[], see Fig. 1(c). In 
the MA model all transitions to neighbor sites have identical rates, if the final site has a lower 
energy than the initial site. Transitions that lead to energetically higher sites require thermal 
activation, 
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The model was introduced for electron or hole transport in amorphous semiconductors. For 
this problem a distance-dependent factor appears in the transition rate. It seems reasonable 
for particle diffusion to restrict the modeling to nearest neighbor transitions. 

For completeness also the model of randomly blocked sites (RBS) is considered ||, cf. 
Fig. 1(d). Particles are allowed to make transitions between accessible sites, with a uniform 
transition rate V . If the concentration p of open sites becomes smaller than the percolation 
concentration p c , long-range diffusion becomes impossible. We restrict the discussion to con- 
centrations p > p c where an asymptotic diffusion coefficient exists. 

3. Previous Results for Single Particles 

What is known about the asymptotic diffusion coefficients of single particles in these different 
models? An exact result is known for the RB model in d = 1, cf. |7H, 
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The curly brackets in Eq. denote the disorder average over the distribution f(T) of the 
transition rates. For simplicity the lattice constant has been set to unity. A consequence of 



Eq. § is that the diffusion coefficient vanishes when the average diverges. This is the case for 
an exponential distribution of barrier heights, v(E) = cr _1 exp(— E/a) at temperatures where 
k B T < a. 

For the RB model in higher dimensions, one has to resort to approximations. One pop- 
ular approach is the effective- medium approximation (EMA). An effective transition rate is 
determined from a self-consistency condition, which reads [|[ f| 



eff 



(3) 



where z is the coordination number of the lattice. Again the average is taken over the distri- 
bution of the rates f(T). The diffusion coefficient is then identical to T e ff. Usually the EMA 



gives good results in d = 3 for not too strong disorder |10 



Another estimate of the asymptotic diffusion coefficient is provided by the critical-path (CP) 



approach of Ambegaokar, Halperin, and Langer ||11|| . In this approach the path from one to 
the other side of a large lattice is considered that overcomes the lowest possible barriers. The 
highest barrier in this path determines the asymptotic diffusion coefficient, and it is related to 
the bond percolation threshold. The CP approach is appropriate for strong disorder in d = 3 



E|. For simple square lattices EMA and CP yield the same result, for a uniform distribution 



of activation energies between and E c , 

D = r "** (4) 

This is in fact an exact result, as a consequence of the self-duality of this lattice ||. 
An exact result is known for the RT model in all dimensions. It reads 
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where hypercubic lattices have been assumed and the lattice constant is unity. This result 
was first derived in ||13|| . That it is exact, can be shown by establishing the strict linearity of 
the mean-square displacement for all (positive) times, for equilibrium initial conditions ||14|| . 
The physical reason for this result is the absence of any time-dependence in the correlations of 
consecutive transitions. Equation 5 follows then easily from the short-time behavior. 

For the MA model to our knowledge no exact result has been published for d = 1. The 
interest in this model concentrated on the physically relevant dimension d = 3. Here EMA 
results have been developed, see the review fT5fl . 

Finally, for the model of randomly blocked sites an approximate result for small concentra- 



tions of blocked sites has been derived by Tahir-Kheli [16]. It reads 



P=r ( 1 - L r+---) ( 6 ) 

where p is the concentration of accessible sites and / the correlation factor for tracer diffusion 
in lattice gases, in the limit of concentration c — > 1. Equation represents an expansion in 



1 — p and the next term is known in d = 2 [17 . Numerical simulations [18] demonstrated that 



Eq. ^ is a very good approximation in d = 3 for p <; 0.8. 
4. Exact Result for d — 1 

Recently an exact expression for the asymptotic diffusion coefficient in one-dimensional 
lattices has been derived fl], [|, which is valid for arbitrary forms of the disordered transition 



rates. The assumption has to made that an equilibrium state exists for the site occupancy in 
the limit of infinite chains. One approach for the derivation of the expression is to consider the 
mean first-passage time (MFPT) of a particle from a starting-site to a terminal site N on a 
segment of a disordered chain Q. An exact expression is known for the MFPT i ^, for fixed 
disorder, in which all transition rates on the segment appear explicitly fT9f . 



N-l 1 N-2 1 N-1 i -n 

'.^E^ + E^E (7) 

k=0 1 k < k + 1 k=0 1 k > k + l i=k+l j=k+l 1 

If an equilibrium state exists, the relation of detailed balance is valid between neighbor sites, 

exp(— f3E{) 

p^ = p^ with * = {eM-vm ( ) 

and (3 = (/c^T) -1 . Here pi is an occupation factor, which is proportional to the occupation 
probability of site i. Equation |] can be used to simplify the second double sum in Eq. [F| One 
obtains as the disorder average of the MFPT in the long segment limit N — > oo 

Comparison of Eq. [| with the inverted asymptotic relation between time and mean-squared 
displacement t = (2D)~ 1 {x 2 } yields the final result 

D = i^r-) ■ (10) 

By analogous derivations as sketched above one can show that the mean-squared deviation 
of the mean-first passage time behaves as 

&} - {W} 2 oc iV 3 . (11) 

This result means that the relative deviations from the mean value become small for large N. 
In this sense it is justified to use the inverted connection between time and distance-squared to 
deduce the diffusion coefficient. 

Equation [10] has been derived in 0, [| from the linear response of the mean particle current 



to an external driving force, invoking the Einstein relation between mobility and diffusivity. 

The physical significance of the result Eq. |10| is that the diffusion coefficient on disordered 
linear chains is determined by transition rates that are weighted by the thermal occupation of 
the initial sites. Loosely speaking, the diffusion coefficient follows from the actual transition 
rates of the particles. 

It will be first verified that Eq. [ll] reproduces previous exact results. This is trivial for the 
RB model, where all pi — 1, and Eq. [| is immediately obtained. Since particles in the RT 
model require thermal activation that is given by the energetic depth of the sites, measured 
from a common origin, the occupation factors are inversely proportional to the transition rates, 
Pi ~ iy . Equation ITU] then reduces to 



D = r {exp(-PE)}-\ (12) 



In view of the equivalence between pi and Tj \ Eq. is completely equivalent to Eq. |5| The 
expression Eq. [12] provides a new interpretation of the general result Eq. [5] for the RT model. 



In this model, the diffusion coefficient is reduced in accordance with the distribution of the 
occupation factors. If these are widely distributed, the diffusivity is reduced accordingly 

A nontrivial application of Eq. |H] is provided by the MA model, in the form presented in 
Eq. |l|. The resulting diffusion coefficient is 



D = -Toiexpi-PE^iexp^Ef)}-^. 



(13) 



The second average is actually a double average over distributions viEi) and viEf). Fig. 2 
gives a numerical verification for the discrete two-level distribution of transition rates. 
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Figure 2: Diffusion coefficient of single particles in a two- level MA model in d = 1 with c± the 
concentration of lower sites. The up transition rates are 0.1 (O), 0.05 (+), and 0.01 (□) in 
units of r . 



5. Effective-Medium Approximation for d > 1 

To derive diffusion coefficients of particles in higher- dimensional disordered lattices one has 
to resort to approximations. Equation [10] gives a suggestion how to handle different forms of 
disorder, when thermal occupancy of the sites plays a role. The suggestion is to use transition 
rates that are weighted by the thermal occupation factors. In the further derivations the 
following rates will be used 

IT = PFir (14) 

The transition rates Eq. |14] are labeled "sym" because they are symmetric in view of the 
relation of detailed balance. This is of direct relevance, since the self-consistency condition Eq. 
|3] requires symmetric transition rates. The use of Eq. [H] in the self-consistency condition Eq. 
|H constitutes the EMA for the different models of disorder. 

The RB model is trivial in the sense that all Pi = I and the self-consistency condition can 
be evaluated with the original symmetric rates. Various results for the RB model have been 
published in the past, see, for instance [f2(J. Evaluation of the self-consistency condition Eq. |3] 
for the RT model reproduces the dimension-independent exact result Eq. |[ A nontrivial ex- 
ample is provided by the MA model, where results for a Gaussian distribution of site energies 



in d = 3 are presented in Fig. 3. Deviations between the EMA and the numerical simulations 
are visible for larger disorder. 
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Figure 3: Diffusion coefficient of single particles in the MA model in d = 3 for Gaussian 
distributions of the energies. The width a of the Gaussian is given in units of fcgT. 

6. Combination of random barriers and random site energies 

In various amorphous materials the observed diffusion coefficients exhibit Arrhenius behav- 
ior, i. e. In/} vs 1/T shows linear behavior. This is surprising because the different models 
discussed above yield curvature in an Arrhenius plot of D. The RT model leads generally to 
downward curvature, an example will be given below. The RB model leads to upward curvature 
in d = 3, an example is given in P0 |. Limoge and Bocquet ]2T] suggested that in a combination 
of the RB and the RT model cancelation of both effects may occur, resulting in an Arrhenius 
behavior of the diffusion coefficient. The methods outlined above are well suited to treat this 
problem. 

Equation |l^ shows that the combination PiT^ enters into the theory. The energy to be 
surmounted from a site over a barrier leading to the neighbor site is the sum of a site energy 
and a barrier energy. The two random quantities are assumed to be independent. In the 
combination pT, the contributions of the site energies cancel each other, as for the pure RT 
model, while the barrier energy remains unaffected. As a result, the diffusion coefficient of the 
combined model is ^ 

D com b = ttD rt D rb . (15) 
1 o 



The coefficient Drt is given by Eq. [12| in any dimension. Drb can be evaluated exactly in 
d = 1, cf. Eq. |2|, and approximately from the EMA in higher dimensions. 

Since the expressions Eqs. |2| and [| are identical in d = 1, the effect of random barriers is 
the same as that of random site energies, if they have the same distribution. As a result, no 
compensation effect will occur in d = 1, and the downward curvature of D in the Arrhenius 
plot of the RT model will be enhanced by random barriers. Compensatory effects are possible 
in d = 3, if the strength of the disorder is adjusted accordingly. 
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Figure 4: Diffusion coefficient of single particles in the combined RT and RB model, for uniform 
distributions of the energies. Numerical results for the RB (+), RT (□), and the combined 
model (O) are shown together with the EMA results. 



The conclusion can be drawn that cancelation of the effects of random site energies and 
of random barriers is possible, under the assumption of independent distributions of the two 
quantities. However, the cancelation requires a specific choice of the parameters of the distri- 



butions. More details for this model will be published elsewhere [22 



7. Many Particles: Previous Exact Results 

Now many particles in lattices with disordered transition rates will be considered. Multiple 
occupancy of the sites will be forbidden, otherwise the particles do not have interactions. The 
derivations will be restricted to the coefficient of collective diffusion, which appears in Fick's law. 
That is, the decay of density deviations is considered, also in the simulations, on a mesoscopic 
scale. If this decay is exponential for longer times, with a coefficient inversely proportional to 
the wavelength squared, the coefficient of collective diffusion can be deduced. 

There are only a few exact results available for the collective diffusion coefficient of lattice 
gases on disordered lattices. The most prominent result is that the coefficient of collective 
diffusion agrees with the diffusion coefficient of single particles whenever all transition rates 
between two neighbor sites are symmetric. There is a cancelation of the exclusion effects of 
double occupancy, leading to D coU = D s p fp3fl. This result holds for the RB model and the 



model of randomly blocked sites. A numerical verification of it was made for the RBS model 
in |Tj|. A consequence of the result is that the derivation of D co u becomes nontrivial for the 
RT and MA models, and for combinations involving the RT model. A trivial limiting case is 
the one of low concentrations of particles, in which D co u —>■ D s p _. 

Lattice gases in the RT model can be considered, in the limit of concentration c — > 1, as 
models with a low concentration of vacancies in random lattices, whose transition rates into the 
sites are specified. The diffusion coefficient can be calculated exactly in this limit in disordered 
linear chains |24| and the result is 
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(16) 



The result Eq. |16] implies that in a two-level model with "free" and "trap" sites, the diffusion 
coefficient has a minimum for a trap concentration of c t = 2/3, not at the concentration c t — 1 



of sites of low energies. 

Another limit that is well understood is the two-level random trap model for d > 1 when 
the concentration of trap sites is small and the free sites are percolating. If the hopping rate 
out of the trap sites is made very small, and if c > c t , the trap sites are saturated by particles 
that are practically immobile. This is effectively the RBS model, where the collective diffusion 
coefficient is given by the single-particle diffusion coefficient. It was pointed out in Sec. 3 that 
Eq. |6] is a very good approximation for the single particle diffusion coefficient in this model. 



Numerical simulations 25, 26 verified that the coefficient of collective diffusion in the situation 



described above is given by Eq. [| 

8. Effective-Medium Approximation at Finite Particle Concentrations 

Apart from the special situations described above, approximate methods are necessary to 
derive the collective diffusion coefficient for the RT and MA models. For arbitrary dimensions 
one would like to extent the EMA to finite particle concentrations. This requires in a first 
step the introduction of effective or mean-field single-particle transition rates. Single-particle 
transition rates can be introduced in two different forms. One possibility is to factorize the two- 
particle probabilities that occur in the master equation and to linearize around the equilibrium 
solution |3] . The derivation of Eq. [16] was based on this procedure. The second possibility is 
to introduce symmetrized transition rates by [[27jj 

-pSym _ Pt(l ~ PjJ£ij 

{P<(l-P,)}" 

Here Pi are the site occupation probabilities in equilibrium, with < Pj < 1. The rates Eq. [17| 
are symmetric as a consequence of the detailed balance condition for lattice gases. They are 
generalizations of the rates Eq. ^ for single particles; note that the Pj and pi are differently 
normalized. The denonimator in Eq. |17] has been introduced for normalization. 

The result for D co u in d = 1 that follows from Eq. [T7| will be discussed in the next section. 
In higher dimensions, the symmetrized rates can be used in the self-consistency condition Eq. 
[| this constitutes an EMA for lattice gases on disordered lattices. Numerical simulations for 
the RT model were compared with the EMA in Ref . p6|, |28| . There is good agreement between 
the EMA and the simulations for smaller c, and rough qualitative agreement for larger c. Also 
the single-particle limit is obtained correctly. 

Results for the MA model are presented in Fig. 5. The exponential distribution was taken 
for the site energies, and the relevant parameter is the ratio between the width a and the 
temperature. This distribution is interesting because the single-particle diffusion coefficient 
vanishes also in the MA model for k^T < a. Due to the saturation of the low-energy states 
a coefficient of collective diffusion exists for finite concentrations. For /c#T < a and small c a 
power-law dependence holds, 



as has been found for the RT model p9| . The simulations agree with the EMA for small c, 
there are deviations at larger c. 

In the limit of infinite coordination number z — > oo, the EMA yields 

DtT = {if m } ■ (19) 

This expression was previously derived as a phenomenological theory of lattice-gas diffusion on 
disordered lattices |30|, [31]. Further details are given in p8[ . 







0.5 

c 



1 



Figure 5: Collective diffusion coefficient in the MA model in d = 3, with an exponential energy 
distribution. The ratio cr/fcgT was 0.7 (O), 1 (+), 2 (□), and 4 (x). The symbols are simulation 
results and the lines represent the EMA. D sp , for a/ksT = 0.7 is indicated on the ordinate. 



9. Open Problems 

Missing parts of a complete theory of diffusion of many particles in disordered lattices 
are evident. Interaction of particles has been neglected here, except the exclusion of double 
occupancy of sites. Also tracer diffusion has been disregarded. Tagged-particle motion involves 
additional correlations, which are certainly difficult to treat in disordered lattices. Disorder has 
been restricted to local defects (sites or bonds); the inclusion of extended defects would be of 
great interest. 

How far is collective diffusion of site-exclusion lattice gases in disordered lattices understood? 
One can say that it is understood from a practical point of view. The major qualitative effect 
is the saturation of trapping sites. Approximate methods exist to deal with the problem in 
higher dimensions. From a fundamental point of view, however, the problem is not yet solved. 
This becomes apparent in d = 1. From the symmetrized rates Eq. |T7| the following expression 
is obtained in d = 1 |27] : 

Doou^^TsyJ- 1 }' 1 . (20) 



An equivalent result follows from the unsymmetrized mean- field rates p4|.The result shows 
strong discrepancies with the numerical simulations at intermediate particle concentrations. 
The discrepancies can be traced back to the neglect of two-particle correlations when the 
single-particle rates are introduced. Garter and Pitis [^] calculated corrections to Eq. p0| ; 
there remain still discrepancies. The conclusion is that the two-particle correlations are not 
understood in time-dependent situations. There remains the task to develop a correct theory 
for this problem, at least in d = 1. 
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